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Abstract 



We prove the existence of at least two geometrically different periodic 
solution with winding number N for the forced relativistic pendulum. 
The instability of a solution is also proved. The proof is topological 
and based on the version of the Poincare-Birkhoff theorem by Franks. 
Moreover, with some restriction on the parameters, we prove the ex- 
istence of twist dynamics. 



1 Introduction and statement of the main re- 
sults 

We are concerned with the equation of the forced relativistic pendulum 



where a is a positive real constant and / is a T-periodic real function with 
mean value zero over a period. 

Results on the existence of solutions to ([T]) and on their qualitative properties 
have been recently given (by different methods) by various authors. We refer 
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to the works of Brezis-Mawhin [5] where it is proved the existence of a T- 
periodic solution and of Bereanu- Torres [2] who provided the existence of a 
second solution not differing from the previous by a multiple of 27r. From 
now on, multiplicity has to be understood in this way. These works use 
very delicate variational techniques. As for topological methods we refer 
to Bereanu- Jebelean-Mawhin [T] where it is proved the existence of two T- 
periodic solutions using degree arguments and upper and lower solutions. In 
this restriction on the parameters a and T is required. 

Another interesting topic concerning ([T]) is the question of the stability of 
the solutions; in this direction, Chu-Lei-Zhang [5] found stable T-periodic 
solutions for the related case of a relativistic pendulum with variable length 
using KAM theory and Birkhoff normal form. 

Our contribution in the present paper contains two main results (Theorem [l] 
and Theorem |3] below) which we describe in what follows. 

Our first result deals with the existence of T-periodic solutions with winding 
number N, i.e. solutions x{-) such that x{t + T) = x{t) + 2N7r,N G Z for 
all t G M. Obviously, T-periodic solutions correspond to the case N = 0. We 
first show (Proposition [T]) that the condition 



(2) 



2A^7r 



T 



< 1 



is necessary for the existence of T-periodic solutions with winding number N; 
this agrees with physical intuition, which suggests (due to the bound arising 
from the speed of light) that we cannot expect solutions for every value of 
N and T. The necessary condition is also sufficient; in fact we prove 

Theorem 1. For every N E 1i such that |^^| < 1 there exist at least two 
solutions such that x{t + T) = x{t) + 2A^7r, N E 7j for all t E M.. Moreover, 
at least one of them is unstable. 

The crucial fact is the passage from the Lagrangian to the Hamiltonian form 
performed by the Legendre transform. This eliminates the singularity and 
allows us to obtain solutions as fixed points of the Poincare map associated 
to the planar Hamiltonian system corresponding to ([T|. If we consider the 
case = we obtain 

(3) q = Hp, p= -Hg 
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where H = ^yp'^ + 1 — acosq — f{t)q. This transform has been used by 
Mawhin in [16] for other purposes. More precisely, in Section 2 we first give 
a version of the Poincare-Birkhoff theorem (Theorem |4]) which is a variant of 
the resuh by Franks |7j. Then, in Section 3 we show (Lemma |4] and Lemma 
|5| that under ^ the assumptions of Theorem |4] are satisfied. The proof of 
Theorem 4.2 in [7] is developed by elaborated techniques from differential 
geometry, while the proof of our Theorem |4] requires more elementary tech- 
niques based on the application of a result of Le Calvez-Wang ^ (Theorem 
|5|. The crucial notion for our argument is that of exact symplectic map (see 
[H] for an introduction on exact symplectic maps). It is worth mentioning 
that Franks gave in [3| applications of his abstract result to the equation of 
the forced (non-relativistic) pendulum; anyway they need some clarification. 
As for a comparison with this and other results on the (non-relativistic) pen- 
dulum, we refer the interested reader to the end of this introductory section. 
For the proof of the instability result in Theorem [T] we use a theorem by 
Ortega [20] , together with the Poincare-Hopf Theorem [S] . 

Our second result provides a restriction on the parameter a which makes the 
Poincare map a "twist map", whose definition is recalled at the beginning of 
Section 2.2. More precisely in Section 3.2 we prove, by a Sturm comparison 
argument. 

Theorem 2. If a < then the Poincare map associated to system ^ is 
twist. If a = vr^/T^ and f does not vanish identically the same result holds 
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true. Moreover the condition a < ^ is optimal. 

This important fact leads us to more information on the number of solutions 
and their stability properties. To state these results, we adapt a definition 
given in [22] saying that a planar first order system in the variables {q,p) 
is degenerate if there exists a curve (g<j(0),ps(0)) such that the application 
s qs{0) is defined from M onto M, satisfies gs+27r(^) = <ls{t) + 2n and 
Ps+27r(^) = Psit), is bijective in [0, 2n) and continuous and for every s G [0, 27r) 
the point (5's(0),Ps(0)) is the initial condition of a T-periodic solutions with 
winding number N. Then we can state 

Theorem 3. // < a < vr^/T^ either the number of isolated T-periodic 
solutions with winding number N is finite or we are in the degenerate case 
and every degenerate solution is unstable. If a = vr^/T^ and f{t) is not the 
trigonometric function asin(^^t), the same result holds true. 
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The qualitative properties of the T-periodic solutions with winding number 
(whose existence is guaranteed by Theorem [l] and, in the context of twist 
maps, by Ortega's result [21] recalled in Theorem [6]) follow from the ap- 
plication of two abstract facts for planar, exact symplectic and twist maps 
(Corollary |2] and Corollary [s] in Section 2) which follow from Theorem [6] as 
well. At this stage, it is important to recall and develop the "intersection 
property" . 

We end this introductory section by a comparison between our results (and 
their method of proof) and some (somehow analogue) results for the forced 
non-relativistic pendulum equation 

(4) X + asinx = f{t) 

with analogous hypothesis as in equation ([T]). It is well known that it has 
at least two T-periodic solutions. The existence of one solution was proved 
by Hamel ^ and rediscovered independently by Dancer [B] and Willem • 
Then the existence of a second solution has been proved by Mawhin and 
Willem [T7] using mountain pass arguments. It is worth mentioning that 
also Franks |7] gave his contribution using his version of the Poincare-Birkhoff 
theorem. More details on Franks's result will be given at the end of section 
3.1. We stress that our approach and results apply also in the case of equa- 
tion (|4]): it is just a matter of computation, the arguments are completely 
analogous. Therefore our results are an improvement also concerning the 
non-relativistic case. 

The paper is organised as follows. In Section 2 we state and prove the 
abstract results for planar maps that we need in the sequel. More precisely, 
in Section 2.1 we deal with the existence of two fixed points of a planar exact 
symplectic map (Theorem |4]); in Section 2.2 we focus on exact symplectic 
and twist maps, recall the existence of two fixed points (Theorem [6]) and 
give some qualitative results on the number of such solutions (Corollary [2]) 
and of their index (Corollary [3]). In Section 3 we prove, using the results of 
Section 2, Theorem [T| Theorem [2] and Theorem [3j In Section 4 we give some 
more detailed results in the particular case when equation ([T]) is autonomous. 
This is done for a better understanding of the main results. 

I want to thank Professor Rafael Ortega who introduced me to this problem. 
Without his constant and patient supervision this work would have not been 
possible. I am indebted to Professor Anna Capietto whose suggestions have 
been very useful for the final layout. 
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2 Two versions of the Poincare-BirkhofF The- 



orem 

2.1 If the annulus is not invariant 

The classical Poincare-Birklioff theorem gives two fixed points of an area 
preserving homeomorphism of the annulus which satisfies a boundary twist 
condition. In this section we will state a variation of this theorem due to 
Franks [7\ in which it is eliminated the invariance of the annulus despite an 
increasing of regularity. 

In a plane with coordinate {6,r), consider, for < a < 6, the two strips 

y4 = M X [—a, a] and i? = M x [—b, b] so that A G B, and the corresponding 
annuh A = §^ x [-a, a] and B = x [-6, b] so that Ac B. 
We will deal with an embedding S : A ^ B defined as follows: 

and 

fe, = F{e,r) 

\n = G{e,r) 

where F, G : — )■ M are functions of class such that 

F{e + 27r, r) = F{e, r) + 27r, G{e + 2-n, r) = G{e, r). 

By the definition of 5* we are allowed to think 5* as a lift of a map S : A ^ B 
defined on the cylinder. Moreover we will suppose that the map S is isotopic 
to the inclusion i, i.e. there exists a map f : A x [0,1] ^ B such that for 
every t G [0,1], ft{x) = f(t,x) is a diffeomorphism from A onto its image, 
fo = S and /i = i. 

Now consider the standard volume form u = d6 A dr and remember that S 
is area preserving if d6 A dr = dOi A dri. Now, if we consider the 1-form 

a = TidOi — rd9, 

we have that S is symplectic if and only if a is closed. But the strip is 
contractible, so closed and exact forms coincide. Summing up we have that 
S is symplectic if and only if there exists V e C^{A) such that 

dV = ndOi - rde. 
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This equivalence is no longer true in the cylinder, in fact the primitive V 
could be multi-valued. So we have the following 



Definition 1. We say that S is exact symplectic if there exists a function 
V e C^{A) such that 



Now we can state the slightly modified theorem by Franks. He dealt with a 
map defined from an annulus into itself, we will need the case of a map defined 
from a strip into itself. His proof deals with very sophisticated techniques 
of differential geometry. We will reach the requested version translating his 
proof in our concrete case of the cylinder so that it will be understandable 
also for people who do not deal with those abstract tools. 

Theorem 4. Consider a map S : A ^ B which is the lift of an exact 
symplectic embedding S : A ^ B isotopic to the inclusion such that S{A) C 
intB. Suppose that the following boundary twist condition is satisfied: there 
exists r > such that 



Then S has at least two fixed points. 

The strategy of the proof is to extend S to an homeomorphism g of the strip 
B. So we can use the fact that the fixed points of the Poincare-Birkhoff 
theorem are in fact fixed points of the lift. More precisely we will use the 
following result in [Ij 

Theorem 5. Let A be a strip and let A be its corresponding annulus. Con- 
sider a map S : A ^ A which is the lift of a homeomorphism S : A ^ A, 
isotopic by homeomorphisms to the identity and area preserving. Suppose 
that the following boundary twist condition is satisfied: 



dV = rid9i-rd9 and V{9 + 2n,r) = V{9,r) V(e,r). 



F{9,a)-9>r, 9e[0,27r) 
F{9,-a) -9 < -r, 9 e [0,2tt). 



(5) 



F{9,a)-9> 0, 
F{9,-a)-9 < 0, 



9 e [0,27r) 
9 G [0,27r). 



Then S has at least two fixed points. 
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At the beginning we will work in the annuli A and B. We observe that in this 
case is not possible to apply the isotopy extension theorem as stated in [181 
p. 63] because B should be without boundary. Anyway, by the hypothesis 
S{A) C intB it is possible to slightly modify it in order to extend S to gi : 
B ^ B so that gi is isotopic to the identity and restricted to a neighbourhood 
of dB is the identity (this is achieved in [TTl Theorem 1.3, p. 180]). Notice 
that gi does not preserve the area out of A. 

Now choose oq slightly smaller than a and define the following subsets that 
we will use during the proof 

io = M X [-ao, ao] 
Ao = §^ X [-ao,ao]- 

In order to apply theorem [5] and get the result, let us prove the following two 
lemmas: 

Lemma 1. It is possible to alter gi : B ^ B finding a diffeomorphism 
g2 : B ^ B such that 

• g2 is area preserving on B, 

• 92\ao = S, 

• 92\dB = Id. 

Lemma 2. It is possible to alter the lift g2 : B ^ B finding g : B ^ B such 
that 

• g is area preserving on B, 

• g has no fixed points out of A, 

• g satisfy the boundary twist condition 

• 9 = 92 on A. 

So theorem |4] will follow from the application of theorem [5] to g. 

Proof of Lemma^ To prove this lemma we will use Moser's ideas, presented 
in ^19j in a more general framework. Let us break the proof in several steps. 
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step 1. Let S+ = §^ X [ao,6] C S. It results 

^(B+) = / det g[{9,r)d9dr. 
Jb+ 

Let gi{9,r) = (6'i,ri). The 1-form ridOi — rdO is exact symplectic on A so 
its integral over any closed path in A must vanish, in particular we have 

/ rde^ ! rde. 

isix{ao} igi(Six{ao}) 

Moreover, because gi is the identity over x {6} we have 

I rde^ I rde 

Js^x{b} Jgi{S^x{b}) 

SO that 

f rd9= I ndOi. 

JdB+ JdB+ 

Notice that 61 = 9i{e,r) and n = ri(^,r) so that d9i = ^d9 + ^dr. 
Finally, by Green's formula. 



B+ 



d9dr= I det g[{9,r)d9dr 

JB+ 



that implies our claim. 



Step 2. Define Vt{9,r) = 1 - det g[{9,r). Then there exist two 
functions a{9,r) and (3{9,r) 27r-periodic in 9 that vanish on dB and such 

that n-^-^ 



Or de- 
Consider the two functions 



a 



{9,r)^- [n{e,r)-n{r)]de and P{9,r) ^ / Q{p)dp 

Jo J ao 



with Q(r) = ^ J^^ Q{9,r)d9. First of all they are of class because gi is 
of class C^. 

Notice that from Step 1 we have that 

(6) / / Vt{9,r)d9dr = 



ao JO 



and, remembering that gi is the identity on a neighbourhood of dB, it results 

(7) Q{e,b)=0 

Moreover, an exact symplectic map is also area preserving so the determinant 
of the Jacobian is 1. It follows that 

(8) n{9,r) = on x[ao, a]. 

By computation we get Q = ^ — ^ on . The fact that a and (5 are 
27r-periodic with respect to 6 is trivial for /3, while comes from (|6]) for a. 
Moreover, remembering ([T]) we have 

a{e,b) = - [ [n{Q,b)-n{b)]de = Q 

Jo 

and by ^ 



and, by Q, a and /3 vanish on §^ x [ao,a]. We can do the same on 
§^ X [—60, — cto] and find a{6,r) and P{6,r) with the same property. Finally 
we can extend these functions to all B setting a{6,r) = and P{6,r) = on 
Aq and the properties on §^ x [ao, a] guarantee the regularity. 

Step 3. Consider the function, for t G [0,1], Vlt{9,r) = (1 — t) + 
tdetg'iiO.r) and define the vector field 

and the associated differential equation 

^ = Xi(t,e,r), f = X2(t,e,r) 

with solution (j)t = (9f, Rt) passing through {9, r) at time t = 0. The solution 
is unique because Xi and X2 are of class (this justifies the hypothesis of 
5" being C^). We claim that 

fi,(e„/?Odet(^^^) = l, tG[0,l]. 
9 



Remember that a map isotopic to the identity is also orientation preserving, 
while the converse is false in the cylinder (as a counterexample take the 
map {0,r) h- )■ {—6,—r)). Hence we have that detg[ > so that the vector 
field is well defined. Notice that if {9,r) G B the solution does not leave 
B because the boundary circles of B are continua of fixed points: it implies 
that (f)t{B) — B. Using Liouville formula for the linearized equation we have, 
for every t e [0, 1] 

Si.(e„ii,)det(?^) 

= n^iOt, Rt) exp{ j tr{-^ + -^){s, e„ Rs)ds} 
= fi,(e„ R,) exp{ i-^a— + - ^^^)ds} 

where in the last equality we used the properties of Q and the fact that 

dnt/dt = -n. 

So, we have to prove that 

R,) exp{ {-^a— + - = 1- 

Passing to the logarithm and differentiating with respect to t this equality is 
true if and only if 

that is the case remembering the definition of 6^ and Rf. 

Step 4. The function §2 = gi ° 0i satisfies the lemma. 
Indeed, by the previous step 

det g'^ = det(^; o <f>,) det <f>[ = n^{e,, R,) det{^^§^^) = 1 

o[U, r) 

that means that is area preserving. Moreover, by the definition of the vector 
field {Xi, X2) we have (j)i\dB = Id and 0i|ao = Id that imply g2\dB = Id and 
5'2Uo ^ S. □ 
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With the same notation let us conclude with the proof of Lemma |2} This 
part does not involve differential forms so we report the version by Franks. 

Proof of Lemma \^ Let g2 '■ B ^ B he the lift (fixed by the boundary twist 
condition) of g2 that extends f : A ^ B. Now consider 

Mo := sup ci(^2(a;), x) 

where d is the distance in and fix M > Mq. So we have that M is greater 
than the distance that a point in B could be moved by g2- Now consider the 
strip y4+ = M X [ao, flQ + e] C v4 such that by the boundary twist condition 
and continuity we have that for all x G A"*" 

(9) PiUx)) - P{x) > r 

where P{xi^X2) = Xi is the projection on the first component. 
Define h : B ^ B hj 

h{e,r) = {e + Mp{r),r) 

where p{r) is smooth, monotone such that p(r) = for r < ao and p(r) = 1 
for r > ao + e. Notice that h is area preserving, it is the identity for r < ao, 
it is a translation by M if r > ao + e, if x G v4+ then h{x) G and 
pCHx)) > P{x). 
Finally consider 

93 = 92 o h. 

For X G v4+ we have, using ^ 

P{gz{x)) - P{x) = P{~g2(h{x))) - P{x) > r + P(h{x)) - P{x) > r > 

which means that we do not have fixed points in A'^. Moreover, if we take 
X = {6,r) G M X [ao + e,b] then gslx) = g2(h{x)) = g2{6 + M, r) and by 
definition of M that means that we do not have fixed points in M x [ao + e, h] 
and the boundary twist condition is satisfied on M x {6}. 
To conclude we consider A~ = M x [— ao — e, — ao] and define analogously 
h'{9, r) = {6 — Mp{r),r) with similar properties of h. Defining g = g-^ o h' we 
get also the complete boundary twist condition. □ 

Let us conclude with a remark on the stability of such fixed points. Re- 
member that a fixed point p of a one-to-one continuous map S : U G — >■ 
is said to be stable in the sense of Liapunov if for every neighbourhood 
V oi p there exists another neighbourhood W G V such that, for each n > 0, 
S'^iW) is well defined and ^"(VT) C V. We have: 
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Corollary 1. If S is analytic, at least one of the fixed points coming from 
theorem^is unstable. 

Proof. For the special case of dimension two, there exists a relation between 
the stability of a fixed point and its fixed point index. In fact it was proved 
in |20] that if a continuous one-to-one map 5* which is also orientation and 
area preserving has a stable fixed point p then either S = Id in some neigh- 
bourhood of p or there exists a sequence of Jordan curves {r„} converging 
to p such that, for each n, 

r„nF^x(5) = 0, i{s,fn) = i 

where r„ is the bounded component of \ r„. 

The set of fixed points can be described by the equation 

{F{e,r) - + {G{e,r) - rf = 

This is an analytic subset of the plane, indeed is the set of the zeros of an 
analytic function. The local structure of these sets is described in [12] : they 
can contain arcs, isolated points and points with a finite number of branches 
emanating from them. In the cases of a non isolated fixed point we can not 
find a sequence of Jordan curves converging to the point and not crossing 
the set, so we have instability. In the case of isolated fixed points, remember 
that the Euler characteristic of the strip is null, so, using the Poincare-Hopf 
index formula for manifolds with boundary [HI p. 447 Theorem 3.1 and p. 233 
Proposition 4.5], at least one fixed point does not have positive index and so 
it is unstable. 

□ 

2.2 The twist condition 

In this section we will consider the case in which the map S satisfies a twist 
condition. More precisely, using the same notation of the previous section, 
consider a map 5* : — )■ M where Q = {{0,r) G : a < r < tlj{6)}, a is a 
fixed constant and ip : M — j-ja, +oo] is a 27r-periodic, lower semi-continuous 
function. 

Definition 2. The defined above map S satisfies the twist condition if 
dOi /dr > 0. If the map is defined from an annulus into itself, we say that 
satisfies the twist condition if the lift satisfies the same condition. 
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Note that the just introduced twist condition and the boundary twist condi- 
tion used in Theorem |4] and [5] are independent. Indeed, consider the annulus 
A = [0,27r] X [-1, 1]: the map 

'ei = e + e 

n = r 

satisfies the twist condition but not the boundary twist one; the map 

6»i = 6* + + |r 
n = r 

satisfies the boundary twist condition but not the twist one. Moreover the 
two conditions coexist in the map 

'01= e + r 

ri = r. 

The twist condition is important in itself because it is fundamental for the 
application of the Aubry-Mather theory. In our case it will lead us to more 
results on the equation of the forced relativistic pendulum. 
It can be proved that 

Theorem 6. Assume that S is exact symplectic and satisfies the twist con- 
dition. Fix an integer N and assume that for each 6 E M. there exists 
rg G]a, ip{9)[ with 

(10) F{e,a) <e + 2N7r < F{e,re). 

Then the system 

'F{e,r) = e + 2N7r 
0(6, r) = r, 

with 9 G [0, 27r[, (0, r) G Vl, has at least two solutions. 

We will explore some consequences of this theorem but let us first point 
out some remarks on the concept of intersection property. 
Consider the cylinder C = §^ x M with the usual covering map n(^, r) = 
(6',r); by the periodicity of S, we can define a new map, also denoted by 
S*, mapping liiVL) into C. Consider also a non contractible Jordan curve 
r in C with positive orientation. The curve F divides C in two connected 
components and let us call -R_(r) the lower one and -R+(r) the other one; 
similarly the curve Fi = S{T) divides C in two connected components and 
let us call _R_(ri) and R+(Vi). We can give the following 
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Definition 3. Using the notations given above, we say that a map S : C ^ 
C has the intersection property if for every non-contractible Jordan curve 
r C C that is the graph of a function, 

rn5(r) ^ 0. 

A map S : C ^ C has the strong intersection property if it has the intersec- 
tion property and for every non-contractible Jordan curve T (Z C that is the 
graph of a function, either T = S{r) or 

s(r) ni?+(r) ^ and s(r) ni?_(r) ^ 0. 

In the proof of tlieorem |6] it was proved that if a map S : C C 
is exact symplectic and preserves the orientation, then it has the strong 
intersection property. Note that the strong intersection property imphes 
the intersection property, while the converse is false. As a counterexample 
consider the following map of the cylinder in itself: 

\ri = r + il){9) 

where ip{9) is a non-negative continuous function such that exists 9^ such 
that ^{9^) = 0. 
Moreover we have that 

Lemma 3. // a map S : C ^ C has the strong intersection property, then 
for every non-contractible Jordan curve F C C that is the graph of a function 
we have 

#{pGrn5(r)}>2. 

Proof. The case F = S'(r) is trivial. In the other case the strong intersection 
property implies that there exist two points p+ G -R+(r) and p- G -R-(r) 
that are connected by an arc 7 C 5'(r). But 5'(r) is a Jordan curve so there 
must exists another arc 7* C ^(r) connecting p_ to J9+ and crossing F in a 
point p* ^ p. □ 

Now we are ready to prove two corollaries of theorem [6] 

Corollary 2. If in theorem^we require also that S{9,r) is analytic then the 
set of the solutions of system |77p is either finite or the graph of an analytic 
2TT-periodic function. 
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Proof. By (10) and the twist condition, we get that for each 6 the equation 



(12) F{e,r) =9 + 2Nn 

has a unique solution r := (p{9). By the uniqueness we have that is 
27r-periodic. Moreover, because of the twist condition we can apply, for a 
fixed 6, the analytic version of the implicit function theorem and get an 
open neighbourhood Uq and an analytic function : f/^ — i- M such that 
F{e,4){e)) = e + 2NTr on Ue. But, by uniqueness, (f){e) = on Uq. Re- 
peating the argument for each 9, we get that is also analytic. So, S{9, (p{0)) 
is the graph of an analytic function in the cylinder: let us call it <f)i{0). This 
comes from the analyticity and the periodicity of 5* and the fact that (j)[6) 



satisfies equation (12). So, by Lemma 3l and 0i must intersect in at least 



two points of the cylinder that are the solutions of system (11 ) when elevated. 
Moreover, from the theory of analytic functions, we know that either the set 
{9 e [0, 2tt] : (f){e) = 01 (^)} is finite or 0(^) = 0i(^) V^. 

□ 



Corollary 3. Suppose that in theorem^ condition (10) is satisfied for some 
N eTj. Let {6,f) be an isolated solution of system (11) and define the map 
T{6,r) = {6 + 27c,r). Then i{T~^S, {6,r)) is either —1 or or 1. 

Proof. First of all notice that z(T~^5', {9,r)) is well defined because {9,f) is 
an isolated fixed point of T'^S*. To compute the index remember that 

2(T-^5, {e, f)) = deg(T-^5 - Id, Bs{e, f)) 

where deg indicates the Brouwer degree and 6 could be chosen small enough 
by the excision property. So we will deal with the degree of the map 

{T~^S - Id){e, r) = {F{e, r) - 2Nn - 6, G{0, r) - r) 
■.= iFie,r),Gie,r)):=Sie,r) 



and to compute it we will use a technique by Krasnosel'skii [13] that allows 
to reduce the dimension. 

By the hypothesis, the point {9, r) is an isolated zero of F{9, r) and by 
the twist condition = > 0. So we can apply the implicit function 
theorem to the equation F{6, r) = and find a function 0(6') defined on 
a neighbourhood of 6 such that F{9, 4>{6)) = and 0(6*) = f. Hence it is well 
defined the function 
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which has 9 as an isolated zero. Now consider the homotopy 

We claim that it is admissible i.e. f) is an isolated zero for every A. Indeed 
consider the system 



(13) 



rAF(0,r) + (l-A)(r-0(0)) = 
\\G{e,r) + {l-\)^{e) = 0. 



Because of the twist condition, if we define J^{6, r, A) = \F{9, r) + (1 — A)(r — 
4>{9)), we have 

— (^,f,A) = A— (^,f) + (l-A)>0 

and so we can apply the implicit function theorem to solve the first equation 
in a neighbourhood of 6 and by the uniqueness the only solution is r = (j){6). 
Substituting it in the second equation we get, because of the definition of 

XG{e, (j){e)) + (1 - A)$(0) = ^ A$(0) + (1 - x)^e) = o ^ $(0) = o 

that, remember, has 9 as an isolated solution. So {9, r) is an isolated solution 
of system (13) and we can choose S small enough such that {9,f) is the only 
solution in Bs{9,f). So we are led to the computation of the degree of the 
map 

{9,r)^ir-m,m) 

that, if $'(0) 7^ 0, can be easily computed by linearization. 

However it could happen that $'(6') = 0. So consider the other homotopy 

where 9 is fixed. To prove that it is admissible, consider the system 

( XF{9,r) + {1 - X){r - (l){9)) = 
\<!?{9) =0. 
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By the definition of ^{6) we liave tliat 6 is an isolated solution of the second 
equation that, substituted in the first one, gives XF{9, r) + (1 — A) (r — f) =0. 
We have that f is a solution and is also the only one, because by the twist 
condition we have 

^[AF(^>) + (l-A)(r-f)] >0. 

So {9, f) is an isolated solution of the system, the homotopy is admissible 
and we can compute the degree of the function 

To use the factorization property of the degree consider the function L{x,y) = 
{y,x). We have 

degiLoW,Bs{e,r)) = deg{L,Bs{0,0))deg{W,Bs{e,f)) = - deg{W, Bs{e,f)). 

Now, by the factorization property 

deg{W, Bs{9, f)) = - deg(F, /,) deg($, I§) 
OF - 

The function $ is defined in dimension 1 so its degree can be either or 1 or 
— 1. Finally i{T~^ {0,f)) can be either or 1 or —1. □ 

Remark 1. An intuitive idea of when these cases could occur is given by 
figure [I| 
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(e) i = -1, $'(0) ^ 
Figure 1: Possibilities of intersections 
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3 Applications to the Poincare map induced 
by a forced relativistic pendulum 

3.1 Proof of theorem [U 

We are looking for T-periodic solutions with winding number N of the equa- 
tion 



(14) 



d 



X 



dt Vv^ 



+ a sinx = /(t), 



under the hypothesis previously stated in the introduction. 
First of all notice that physical intuition suggests that it should not be pos- 
sible to have such solutions for every and T, because of the bound given 
by the speed of light. This is a necessary condition that holds for a larger 
class of equations, namely: 

Proposition 1. Let x{t) be a T-periodic solution with winding number N of 



(15) 



d 



X 



dt \^/T 



X^ 



F{t,x) 



where F{t,x) is continuous and T-periodic in t. 
Then 



(16) 



2N7T 



T 



< 1. 



Proof. By Lagrange theorem we get 



|2A^7r| 



\x{t + T) -x{t)\ = \x{c)T\ 



for some c G (t, t + T). But the domain of equation (15) is \x{t)\ < 1 for all 
t, so 

|2A^7r| < T. 

□ 



In this section we will see why the relativistic condition (16) is also sufficient 
to have T-periodic solutions with winding number A^. The proof will be an 
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application of theorem [4] considering S the Poincare map. 



First of all, let us perform the change of variables 
(17) 



y{t) = x{t) - —t. 



Notice that in this way y(t + T) = y{t) and T-periodic solutions with winding 



number of ( 14 ) correspond to classical T-periodic solution of 



d 
dt 



y + 



2Ntt 



2Ntt ' 

+ asm\y + -j^t 



m. 



We will find T-periodic solutions of equation (18) as fixed points of the 
Poincare map. 



Equation (18) can be seen as the Euler-Lagrange equation coming from the 
Lagrangian 



L{y,y,t) = -Jl-{y + —^y + acos{y + —^t] + f{t)y 



T 



T 



and presents singularities. On the other hand, if we perform the change of 
variables given by the Legendre Transform 



(19) 



we get the Hamiltonian 



(i = y 

p-dL 



dy 



H{q,p,t) = [py-L{q,x,t)\y=y^p) = Vr + l-^^p-acos{q+^^t)-f{t)q 
and the new Hamiltonian system 

p 2Ntt 



(20) 



q 

.p 



-a sin (q+^t) + fit). 



that has no singularities and has solutions globally defined. 

Now let S{qo,po) = iQiqo,Po), P{qo,Po)) = {q{T,qo,po),p{T,qo,po)) be the 
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Poincare map associated to system (20) that is well defined because of the 
boundedness of the second term. First of all notice that the Poincare map is 
isotopic to the inclusion. The isotopy is simply given by the flow 0, in fact 
we have the isotopy 

0(tT,go,Po) te[0,l]. 

Notice that this isotopy is valid also in the cyhnder because by periodicity 
we have 

0(t,go + 27r,po) = 0(t,go,Po) + (27r,0). 
Moreover, inspired by [13] we can prove 

Lemma 4. The Poincare map S is exact symplectic in the cylinder. 
Proof. Let us call 9 = qo,r = po and K = Consider the function 



V{9,r) 



P 



K)-H{q,p,t)]dt 



l+p^{t,e,r] 



+ a cos(g(t, 9, r) + Kt) + /(t)g(t, 9, r)]dt. 



First of all, it follows from the periodicity of (18) and the change of variables 
(|19|) that q{t, 9 + 27r, r) = q{t, 9, r) + 2tt and p{t, 9 + 2tt, r) = p{t, 9, r). Hence 
we have 

V{9 + 2TT,r) = V{9,r) 

using the hypothesis of the null mean value of /. Now let us compute the 
differential dV. We have 



(21) 



P 



dv da 



p 



dp .dq. 
+ P7^]dt 







using the second equation in (20). Now, integrating by parts and using the 

dp 



first equation in (20) we get 



/ At - b^f - I v%t - b^f 



P 



{l+p^)y^89 
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that, substituted in (21) gives 



Analogously we can get 

K=P(T)|(T)-M0)|(0). 

Hence dV = pidqi — pdq, that means that the function V will satisfy the 
thesis. □ 

Finally, the relativistic effect will give also the boundary twist condition: 

Lemma 5. // |^^| < 1 then there exists p > and r > such that 

Q{q, ~p) — q < —r and Q{q,p) — q > r. 

Proof. Let us prove the first inequality, being the second similar. Let us call 

Ntt 
T ' 

p 



K := < 1 and consider the function, coming from system (20), 



A{p) 



We have that A{p) is an odd increasing function such that ^(0) = and 
\mip^±oo A.{p) = ±1. Since li^l < 1, by continuity, we can find p > such 
that 

j A{p) > K ioT p > p 
\A{p) < K ioi p < —p. 

Now, integrating the second equation of (20) we get, for t G [0,T] 

p{t) = Pq — / asiii{q{s) + K)ds + I f{s)ds<po + t{a + 
Jo Jo 

and so we can find p > large enough so that if po < —p then p{t) < —p for 
t G [0, T]. It means 

pit) 2A^7r , 

that is q{t) is decreasing if t G [0,T] so, 

QiQo, -p) = q{T, go, -p) < g(0, go, -p) = go- 
Now a standard compactness argument concludes the proof. □ 
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Now it is straightforward the apphcation of theorem |4] choosing the strip 



y4 = M X [— p, and the fact that solutions of system (20) are globally 



defined implies that we can find a larger strip B such that S{A) C intB. 



Finally we have that the right-hand side of (20) is analytic in then, by 
analytic dependence on initial conditions, also the Poincare map is analytic. 
Notice that we do not need the analyticity of / [111 p. 44]. So, using corollary 
[T]we get the instability of one solution. 
Then theorem [T] is proved. 

Remark 2. Similar results on the classical pendulum have been obtained 
by Franks in |^ Proposition 5.1]. He proved the existence of fixed points 
for the Poincare map using his version of the Poincare-BirkhofJ theorem 
and affirmed that they should have positive or negative index. This result 
needs some clarification. In fact there is another possibility: there could 
be only a continuum of fixed points and the fixed point index could not be 
defined. Consider the equation of the classical pendulum: the existence or 
not of forcing terms f of null mean value such that the periodic solutions 
are represented only by a continuum in still an open problem. Anyway, as a 
related example consider the equation 

27r , 

y + asm{y + —t) = 

where the potential depends on time. Its T -periodic solutions correspond, via 
the change of variables x = y + ^t to solutions x(t) of 

x + a sin x = 

such that x{t + T) = x{t) + 2tt. These solutions forms the graph of a function 
in the phase space, so it is impossible to define the index. 

3.2 Proof of theorem [2] and consequences 



We will prove the result for the Poincare map of system (20). In the case 
a = TT^/T^ we will have to prevent the function / from being the trigono- 
metric function asin(^^t). 

According to the previous section let us call S{qo,po) = {Q{qo,Po), P{qo,Po)) = 



{q{T, qQ,pQ),p{T, qo,Po)) the Poincare map associated to system (20) in which 
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we suppose a < vr^/T^. We have to prove that -^{T, qo,Po) > 0. 



If we call 



dg dp 

- {t, go, Po) y{t) = Po) 



we know from the elementary theory of ODEs that the vector {x{t),y{t)) 
satisfies the variational equation 

1 




(l+p2(t,go,Po))3/2 



y 



-acos(g(t,go,Po) + ^^t)x 

1 



that is equivalent to the problem 

d 



(22) 



^^(a;(l +p2(t,go,po))3/2) + acos(g(t,go,Po) + ^t)x 



x{0) 
x{0) 







i_\3/2 



Now consider the equation 
(23) 



71 



and first suppose that a < vr^/T^. In this case we have that 



2N7T 



(1 + p{tYf/^ > 1 and a cos(g + -^t) <a< n^/T 



T 



then (23) is a strict Sturm majorant of (22). So the Sturm theory and the 
fact that the function z{t) = sin(t^) is a solution of (23), prove that x{T) > 
and the thesis will follow. 

Now consider the case a = vr^/T^. First of all the hypothesis f{t) ^ 
asm(^j^t) prevents q = 2kTi from being a solution. This means that there 
exists an open subset of positive measure of [0,T] on which q ^ 2k7T and so 



2^2 



COS q(t)dt < 



T ^ 



dt. 



In this case we can use a generalization of the classical Sturm separation 
theorem. It can be achieved adapting the classical proof (cf. [^I^j) to our 
framework. Consider the argumentum 9i and 62 respectively of (22) and 
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(23) coming from the Prufer change of variables; then we can conclude 
that 0i{T) > 02{T). Remembering that in this framework we have that 
x{t) = 6{t) = for some G Z and we are rotating in the clockwise 
sense, we can conclude using the same argumentation of the previous case 
translated into the phase-space {x,px). The proof of theorem [2] is complete. 



2 2 

Remark 3. The condition a < ^ is optimal. Indeed suppose a > ^ and 
consider the autonomous system 

p= -asm{q + ^t). 

Notice that {p = 0,q = —^^t) is an obvious solution. As before consider the 
variational equation 

^ix{l + p'^{t, qo,po)f^'^) + acos{q{t,qo,po) + = 0. 

at I 

Notice that evaluated in the above solution it is nothing but 

X + ax = 0. 

2 

Using Sturm comparison with y + ^y = we can conclude analogously as 
before that x(T) < 0: it means that we do not have the twist condition. 
Finally note that in the case a = we have x{T) = and again the twist 
condition fails. 

The fact that the Poincare map satisfies the twist condition allows us to 
enter in the huge chapter of twist maps. In particular we will get some more 
results on equation ([T]). To state it remember that a T-periodic solution with 
winding number is said to be isolated if there exists 6 > such that every 
solution {q(t),p(t)) satisfying 

0< |g(0)-g(0)| + |p(0)-p(0)| <5 

is not T-periodic with winding number N. 
We have 

Theorem 7. // < a < vr^/T^ either the number of isolated T-periodic 
solutions with winding number N is finite or we are in the degenerate case 
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and every degenerate solution is unstable. Moreover, in the first situation, 
the index of such solution is either —1 or or 1. If we consider the case 
a = TT^/T^ and we add the hypothesis that f{t) is not the trigonometric 
function asin(^^t), we get the same results. 



Proof. By the twist condition we can apply the results of section |2.2[ In 
particular theorem ^ runs with = {(g,p) G : —p < p < p} where p 



comes from Lemma [5] Indeed, if we take rg = p — e with e sufficiently small, 



condition (10) holds with = by continuous dependence, and from the 
previous section we have that the Poincare map is exact symplectic. This is 
another way to find two periodic solutions. Notice that it is a weaker result 
because we need the restriction on the parameter a. 

Anyway the Poincare map is analytic, so, by corollary [2| we have that fixed 
points either are isolated or form the graph of an analytic 27r-periodic func- 
tion. Moreover by corollary |3] we have the informations on the degree. 
The translation of these results from the Poincare map to the differential 



equation gives informations on the periodic solutions of system (20) and 



by the change of variables (17) we get analogous results on the T-periodic 



solutions with winding number of system ( 14 ) . □ 



4 The autonomous case 

Finally, consider the case / = 0, i.e. the autonomous equation 

(24) ^ ( , ^ H +asinx = 

2 

with a < that can be treated with a phase portrait analysis. Let us 
consider the case T = 2n, so that a < 1/4. 

First of all it is easily seen that the points {kir, 0), A; G Z are constant solutions 
in the phase space {x,x). This analysis is quite simple because the energy 

(25) E(x, x) = 1^ _ ^ — a cos x + a 



X 



is a first integral and we suddenly reach the conditions 
(26) E>1 and - 1 < i; < 1. 
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Remembering that a < 1/4 we get the phase portrait in figure |4] where we 
have the constant solution (0, 0) for E = 1, periodic orbits for 1 < < l+2a, 
the heterochnic orbits for E = 1 + 2a and the unbounded solutions for 



E > 1 + 2a. Moreover, from the first integral (25) we can see the velocity as 



a function of the time and energy and a; — )■ ±1 as — ?■ +oo depending on 
the sign of x{0). 




Figure 2: Phase portrait 

Now we turn to the study of the period of the periodic orbits, in particular, 
in order to complete theorem[l]we will study the number of 27r-periodic orbits. 



Proposition 2. The only 2iT-periodic solutions of equation (24) are the con- 
stant ones (fcvr, 0) with /c G Z. 

Proof. We will prove that the period of every orbit (except for the constants 
one) is strictly greater than 27r. To do so, by the symmetries of the phase 
portrait, it is enough to prove that for every non-constant periodic orbit, 
x{n, xq, 0) > with xq < 0. 

Let us write the solution x{t,E) such that x{0,E) = arccos(l — E/a) and 
x{0,E) 

dx 



and compute ff (tt, E). Remembering (26) and a < 1/4 we have 



dE 



> for E >1. 



Notice that the point (0,0) is a strict minimum of E{xq,Xo) an so 
x(7r, E{0, 0)) < x{7r, E{xo, xq)) V(xo, xq) ^ (0, 0). 
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Now, remembering that E is constant on the solutions, we have that for every 
initial condition {xq,xo) such that 1 < E{xo,xo) < 1 + 2a there exists x < 
such that E{xo, xo) = E{x, 0) and x{7i, E{x, 0)) > x^tc, E{0, 0)) = 0. 

□ 

Looking for T-periodic solutions with winding number we can do the 
following. By the phase portrait analysis we got that for E > 1 + 2a the 
solution is unbounded and the orbit in the phase plane is the graph of a 
function. In this case we will show 

Proposition 3. Fix \N\ > 1 and T > 2n such that ^ < 1. Then there 
exists exactly one value of the energy E > 1 + 2a such that 

x{T + t,E) = x{t,E) + 2N'K. 

Proof. Let us prove in the case t = 0, x{0) = and > 0. 
Remembering the energy ( [25] ) we can define a function T^^E) such that 
x(T/v(-E')) = 2Nti (i.e. T]\}{E) is the time needed by a solution starting from 
at t = to reach 2A^7r), namely 

'•^^^ dx 

Tm{E) ■ 



Jl - 



{E+a cos x—a)^ 

Notice that it is continuous, monotone decreasing in E and 
lim Tn(E) = +00, lim Tn(E) = 2N'k. 

The proposition is proved if we can find E > 1 + 2a such that Tn{E) = T. 
It is automatic using the properties just mentioned and the fact that by 
hypothesis T > 2N7c. 

□ 
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